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ENVIRONMENT AND MASS DEPENDENCIES OF GALACTIC A SPIN 
PARAMETER: COSMOLOGICAL SIMULATIONS AND SDSS GALAXIES 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to translate your abstract 
into Spanish, the editors will do it for you. We use a sample of galaxies from the Sloan Digital Sky 
Survey (SDSS) to search for correlations between the A spin parameter and the environment and mass of 
galaxies. In order to calculate the total value of A for each observed galaxy, we employed a simple model of 
the dynamical structure of the galaxies which allows a rough estimate of the value of A using only readily 
obtainable obscrvablcs from the luminous galaxies. Use of a large volume limited sample (upwards of 11,000) 
allows reliable inferences of mean values and dispersions of A distributions. We find, in agreement with some 
N-body cosmological simulations, no significant dependence of A on the environmental density of the galaxies. 
For the case of mass, our results show a marked correlation with A, in the sense that low mass galaxies present 
both higher mean values of A and associated dispersions, than high mass galaxies. This last direct empirical 
result, at odds with expectations from N-body cosmological simulations, provides interesting constrain on the 
mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models. 

ABSTRACT 

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the 
A spin parameter and the environment and mass of galaxies. In order to calculate the total value of A for 
each observed galaxy, we employed a simple model of the dynamical structure of the galaxies which allows a 
rough estimate of the value of A using only readily obtainable observablcs from the luminous galaxies. Use of 
a large volume limited sample (upwards of 11,000) allows reliable inferences of mean values and dispersions of 
A distributions. We find, in agreement with some N-body cosmological simulations, no significant dependence 
of A on the environmental density of the galaxies. For the case of mass, our results show a marked correlation 
with A, in the sense that low mass galaxies present both higher mean values of A and associated dispersions, 
than high mass galaxies. This last direct empirical result, at odds with expectations from N-body cosmological 
simulations, provides interesting constrain on the mechanisms of galaxy formation and acquisition of angular 
momentum, a valuable test for cosmological models. 
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1. INTRODUCTION 

One of the most studied parameters in numerical 
simulations of formation and evolution of galaxies, is 
the A spin parameter, first introduced by Peebles in 
1969, where it was used to test the gravitational in- 
stability picture, as a theory for the origin of galactic 
angular momentum. Since then, there have been sev- 
eral studies using the A parameter as an indispens- 
able tool of analysis, or characterizing its distribu- 
tion in cosmological N-body simulations. In general, 
the distribution of this parameter, coming from nu- 
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merical simulations, is well fitted by a log- normal dis- 
tribution, characterized by two parameters; Ao, the 
most probable value, and a\, which accounts for the 
spread of the distribution. In a recent work, Shaw 
et al. (2006) showed a compilation of results for es- 
timates of these parameters, arising from numerical 
studies performed by several authors, lying in the 
ranges: 0.03 < A < 0.05 and 0.48 < a x < 0.66, 
and some of us in a previous work (Hernandez et al. 
2007), using the same sample of 11 597 galaxies from 
the Sloan Digital Sky Survey (SDSS) we treat here, 
obtained the distribution of this parameter, with re- 
sults well fitted by a log-normal function with pa- 
rameters A = 0.04 ± 0.005 and cr A = 0.51 ± 0.05, 
for the first time derived from a statistical sample of 
real galaxies. 
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With the recent explosion of high resolution cos- 
mological simulations, studies concerning the spin 
parameter begun to explore dependencies of the A 
parameter with other physical parameters, mainly 
the total mass of the galaxy and the environment in 
which the galaxy is immersed. The results are as 
varied as the techniques employed to perform and 
analyse the simulations, but there is good agree- 
ment that for the case of a A-environment relation, it 
should be very weak or inexistent (Bett et al. 2007, 
Maccio et al 2007). The case of the A-mass relation 
presents more discrepancies between results coming 
from different groups, some claiming the existence 
of a weak relation (Jang-Codell & Hernquist 2001, 
Bett et al. 2007) and others not finding any (Maccio 
et al. 2007). The aim of the present work is to de- 
rive distributions of the spin parameter as functions 
of mass and environment, using a first order esti- 
mate of galactic halo A applied to a large sample of 
galaxies from the Sloan Digital Sky Survey (SDSS), 
in order to constrain and infrom theories of angular 
momentum acquisition and galaxy formation. 

2. THEORETICAL FRAMEWORK 

The angular momentum is commonly character- 
ized by the dimcnsionlcss spin parameter: 



A = 



L I E |Va 



(1) 



GM 5 / 2 

where E, M and L are the total energy, mass 
and angular momentum of the configuration, respec- 
tively. In Hernandez & Cervantes-Sodi (2006), we 
derived a simple estimate of A for disk galaxies, us- 
ing a simple model employing two components to 
describe the galaxy, a disk for the barionic compo- 
nent with an exponential surface mass density E(r); 



E(r) = S e 



-r/R d 



(2) 



where r is a radial coordinate and Eo an d Rd are 
two constants which are allowed to vary from galaxy 
to galaxy, and a dark matter halo having an isother- 
mal density profile p(r), responsible for establishing 
a rigorously flat rotation curve Vd throughout the 
entire galaxy; 
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We assume: (1) that the specific angular momen- 
tum of the disc and halo are equal e.g. Fall & Efs- 
tathiou (1980); (2) the total energy is dominated by 
that of the halo which is a virialized gravitational 
structure; (3) the disc mass is a constant fraction of 



the halo mass F = Md/Mn- Introducing a Tully- 
Fishcr relation and fixing constants with Galactic 
parameters, we obtain: 



A = 21. 



Rd/kpc 



{Vd/kms^f/ 2 



(4) 



The above equation allows a direct estimate of 
A for any disk galaxy. The corresponding equation 
for elliptical galaxies was derived in Hernandez et 
al. (2007) modeling two components; a barionic 
matter distribution with a Hernquist density pro- 
file and a dark matter halo showing no difference 
from that of disk galaxies (Treu et al. 2006). Con- 
sidering the structure of the galaxy as virialized, the 
mass is obtained in terms of the half light radius and 
the velocity dispresion a. By dimensional analysis, 
the specific angular momentum is characterized by 
(GMba) 1 / 2 , where Mb is the baryonic mass and a is 
the mayor axis. Introducing a numerical factor to ac- 
count for the dissipation of the baryonic component, 
its dynamical pressure support and the projection 
effects, we finally obtain 



A = 



0.1173e(a/fcpc) 2 / 7 
(cr/fcms -1 ) 3 / 7 



(5) 



The accuracy of equation (4), was tested in 
Cervantes-Sodi et al. (2008), making use of numer- 
ical simulations, showing an agreement between the 
estimate and the real value of A always better than 
30 per cent, in most cases much better. 

3. COMPARISONS OF A DISTRIBUTIONS 
FROM THE SDSS AND COSMOLOGICAL 
N-BODY SIMULATIONS 

3.1. Observational sample 

The sample of real galaxies employed for the 
study comes from the SDSS Data Release 5 
(Adclman-McCarthy et al., 2007). It is a volume 
limited sample having galaxies in the redshifts in- 
terval 0.025 < z < 0.055 and absolute magnitudes 
M r — hlogh < —18.5. Since most of the studies con- 
cerning spin distributions from simulations presents 
their results at z = 0, we limited the sample to 
low redshifts. This sample contains 32 550 galax- 
ies for which Choi et al. (2006) have determined the 
exponential disc scales, absolute magnitudes, veloc- 
ity dispersions, de Vacouleurs radii and seeing cor- 
rected isophotal ellipticities for each galaxy, assum- 
ing a ACDM universe with fl M = 0.27, Q A = 0.73 
and h = 0.71. To obtain A for each galaxy, we need 
to discriminate between elliptical and disc galaxies; 
in order to do that we used the prescription of Park 
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& Choi (2005), in which early (ellipticals and lentic- 
ulars) and late (spirals and irregulars) types are seg- 
regated inau-r colour versus g — i colour gradient 
space and in the concentration index space. In order 
to apply equation 0] to the disc galaxies of the sam- 
ple, we need the rotation velocity, which is inferred 
from the absolute magnitude using a TF relation, 
then, to avoid the problem of internal absorption in 
edge-on galaxies, we employed only spiral galaxies 
having axis ratios > 0.6 and inferred rotation ve- 
locities in the range of 80 < Vr < 430, well within 
the range of applicability of the TF rcaltion we are 
using. After applying these two cuts, and removing 
a randomly selected fraction of ellipticals, so as to 
maintain the original early to late type fraction, we 
are left with a total of 11 597 galaxies. 

For the elliptical galaxies we apply equation [5j 
as fully described in Hernandez et al. (2007), using 
the corrected isophotal ellipticities to calculate the 
eccentricity. 

3.2. Simulation & Subhalo Finding 
To compare our observational sample with a nu- 
merical simulation, we have made a cosmological 
N-body simulation of a ACDM model, Dutton, A. 
A., van den Bosch, F. C, Morre, B., Potter, D., 
& Stadel, J in a cubic box with the side length of 
614/i _1 Mpc. The model parameters are h = 0.7, 
n m = 0.27, Q b = 0.046, and Ct A = 0.73, given by the 
WMAP 1-year data (Spergel et al. 2003). The sim- 
ulation gravitationally evolved 1024 3 CDM particles 
from redshift z = 48 to z = taking 1880 global 
time steps. The linear density fluctuations of matter 
at the present epoch is normalized by as = 0.9, the 
RMS density fluctuation at 8/i _1 Mpc for a tophat 
window function. The mass of each simulation par- 
ticle is M p = 1.4 x 10 10 /i~ 1 M Q and the force spatial 
resolution is about 60h~ 1 kpc. 

To identify subhalos in the simulation, we first 
extracted the dark matter particles located in viri- 
alized regions by applying the standard Friend-of- 
Friend (FoF) method with a linking length equal to 
one fifth the mean particle separation. This length 
is a characteristic scale for the identification of viri- 
alized structures. To each FoF particle group, we 
applied the PSB method (Kim & Park 2006) to find 
subhalos. 

For comparisons with our observational results 
we randomly selected 100 000 halos from the simula- 
tion, for which the spin parameter A was calculated 
numerically. The normalized local density parame- 
ter is obtained following the same method as used 
for our SDSS sample (for more details see Park et 
al. 2007). 



4. RESULTS 

To study if there is any correlation between spin 
and the environment, we use an estimation of the 
local background density defined directly from ob- 
servations, that is continuous and able to character- 
ize the full range of galaxy environments. It mea- 
sures the local number density of galaxies at a given 
smoothing scale (Park et al. 2007). We divided both 
of our samples, the SDSS and the simulation, into 6 
cuts according to the normalized environmental den- 
sity pj < p >, and for each cut we obtained the mean 
A value and its dispersion. The result is presented in 
figure 1 left panels, where the value of total infered 
A versus pj < p > for each galaxy is plotted and 
we show the mean A value for each cut as a function 
of the normalized density, the solid error bars corre- 
sponds to the dispersion and the broken error bars 
to the standard error in the calculation of the mean 
value, this convention will be kept in the following 
figure. As can be seen, there is no clear trend, tak- 
ing into account the dispersion in each bin, in neither 
the SDSS nor the simulation sample. 

On the right panels of figure 1 are presented the 
plots of A versus mass for the SDSS sample (top 
panel), and the simulation sample (bottom panel). 
We notice that for the SDSS sample, the mean value 
of A tends to decrease as the mass increases. But 
not only the mean value, also the dispersion of A 
tends to decrease as the value of the mass increases. 
The effect is stronger using only the disc galaxies, 
where our estimate is more precise (see Cervantes- 
Sodi et al. 2008). Large galaxies are seen to form a 
more coherent low A sample with small dispersion, 
while small galaxies show mean values of A of about 
a factor of 3 larger than what is found for our most 
massive galaxies. The dispersion about mean val- 
ues shows the same trend, with small galaxies in our 
SDSS sample having a much larger dispersion than 
the large galaxies. The trend is not reproduced by 
the simulation used in this work, and neither has it 
ever been reported, to the large extent seen in the 
SDSS sample, in any cosmological N-body simula- 
tion found in the literature. However, precisely such 
a trend was recently confirmed by Berta et al. (2008) 
using a sample of 52 000 galaxies from the SDSS em- 
ploying a similar method as the described here to es- 
timate A. We take our results as evidence of a mech- 
anism of acquisition of angular momentum for galax- 
ies, where it is the ambient tidal field what torques 
up a halo, with little participation of repeated merg- 
ers. An extensive development of these ideas can be 
found in Cervantes-Sodi et al. (2008). 
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Fig. 1. Relations between galactic spin and environmental density and mass, showing the mean A value for each bin 
as a function of the normalized density, the solid error bars correspond to the dispersion and the broken error bars 
to the standard error in the calculation of the mean value. Top panels corresponding to the SDSS sample using the 
11 597 galaxies, on the left-hand A as a function of environmental density, on the right-hand A as a function of mass. 
Bottom panels corresponding to the sample of 100 000 simulated dark matter haloes, on the left-hand A as a function 
of environmental density and on the right-hand as a function of mass. 
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